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Abstract 

. Real collective density variables C(fc) [c.f. Eq.1.3)] in many-particle systems arise from non-linear 

transformations of particle positions, and determine the structure factor S{k), where k denotes 
^ I the wave vector. Our objective is to prescribe C{k) and then to find many-particle configurations 

' that correspond to such a target C(k) using a numerical optimization technique. Numerical re- 

^ , suits reported here extend earlier one- and two-dimensional studies to include three dimensions. In 

addition, they demonstrate the capacity to control S{k) in the neighborhood of |fc| = 0. The opti- 

■ mization method employed generates multi-particle configurations for which S{k) oc \k\", \k\ < K, 
G\ ; 

^2 . and a = 1, 2, 4, 6, 8, and 10. The case a = 1 is relevant for the Harrison-Zeldovich model of 

I the early universe, for superfluid ^He, and for jammed amorphous sphere packings. The analysis 

' also provides specific examples of interaction potentials whose classical ground state are configu- 

■ rationally degenerate and disordered. 
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1. INTRODUCTION 



Spatial arrangements of particles in many-body systems exhibit wide diversity that arises 
from the interactions that are present, and from the prior history of those systems. One of the 
available analytical tools that has proved useful for describing those spatial arrangements, 
whether for individual cases or for ensemble averages, is the set of collective density variables. 
These are conventionally defined for N identical particles in the following way: 

N 

p{k) =^eKp{ik-rj) (1.1) 

where rj denotes the location of particle j, and the k are the wave vectors appropriate for 
the containing volume and boundary conditions. 

Collective density variables have played an important role in a variety of problems in 
condensed matter physics. Specifically, they have been used in illustrating that large-scale 
density variations in superfluid "'He arc in fact long- wavelength phonons [1]. Introduction of 
the appropriate collective coordinates is natural in describing independent plasma oscilla- 
tions brought about by the long-range Coulomb interactions between electrons in metals [2] . 
Furthermore, application of these variables has aided in the derivation of self-consistent inte- 
gral equations for pair correlation functions in classical fluids [3] , in obtaining corrections to 
the random phase approximation for the electron gas [4] , and in generating classical ground 
states for particle systems [5-7]. But in spite of the fact that collective density variables 
appear widely in the literature [1-10], the nonlinearity of the transformation in Eq. (1.1) 
from particle positions to collective variables entails nontrivial mathematical properties that 
are still incompletely understood [5, 6]. The present paper is devoted to clarifying some of 
the remaining issues. 

Many physical applications of collective density variables focus on the structure factor 
S{k) for the many-particle system involved: 



^2 (1.2) 
= 1 + j^CXk) 



where the real quantities C{k) are the following: 
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AT-l N 

^(^) = E E cos[fe-(r,-n)]. (1.3) 

j=i i=j+i 

In view of the fact that the phase angles of the p{k) are irrelevant for most applications, it 
suffices to focus attention on the C(fc)'s. 

A considerable challenge involves determining what sets of C{k) values correspond to 
attainable particle configurations, and how to generate and describe those special configu- 
rations, including ground-state structures. In particular, it is important to understand the 
extent to which these real collective variables at small wave vectors k are controllable. Al- 
though we begin by considering the general situation, much of the attention in the following 
will involve the examination of "hyperuniform" systems [11, 12], namely those in the infinite 
system limit for which: 

hmS{k)=0. (1.4) 
fe— »o 

This defining characteristic of hyperuniformity states that the usual mean-square particle- 
number fiuctuations increases less than R^, where R denotes the linear size of an observation 
window and d is the space dimension [11]. Considering the fact that various hyperuniform 
physical systems exhibit characteristic k dependence of their structure factors near the 
origin {e.g., the ground state of liquid ^He [13-15] as well as random, jammed hard-sphere 
packings [16] and the early Universe [17]), it becomes important to understand what N- 
particle configurational imphcations stem from these specific S{k) forms. This has guided 
the analysis detailed below. 

Because our objective is to prescribe C{k) [or, equivalently S{k)] and then to find many- 
particle configurations that may correspond to such a target structure function, this problem 
can be regarded to be an inverse problem. The analogous inverse problem in real space in 
which the pair correlation function is prescribed has received considerable attention in the 
last several years [18-22]. This class of inverse problems has come to be called construction 
problems [18-21]. A priori, a prescribed pair structural function is not necessarily realizable 
by a many-particle configuration. A solution to the construction problem therefore provides 
numerical evidence for realizabilty of a target pair structural function. Thus, the present 
investigation has important implications for the realizability problem of statistical mechanics, 
which seeks to determine the necessary conditions that a prescribed pair correlation [or its 
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equivalent Fourier representation, S'(k)] must possess in order for it to correspond to a 
many-particle system [19, 22, 23]. 

The present paper extends our earlier one-dimenisonal [5] and two-dimensional [6] studies 
of collective coordinate control of density distributions. Our focus in these two studies was 
to consider continuous and bounded pair potentials and to constrain the corresponding 
collective parameters C{k), with wave vector k magnitudes at or below a chosen cutoff, 
to their absolute minimum values. In other words, density fluctuations for those fc's were 
completely suppressed. In our two-dimensional investigation [6], we were able to distinguish 
between three different ground-state structural regimes as the number of constrained wave 
vectors were increased - disordered, wavy crystalline, and crystalline regimes. Evidence for 
a disordered or irregular ground state, a counterintuitive notion, had heretofore not been 
provided. In the present work, we extend these results not only to three-dimensional ground 
state problems but to more general two- and three-dimensional hyperuniform many-particle 
systems. 

The next section provides some of the mathematical structure needed to understand how 
fixing values of sets of the real collective variables C{k) exerts control over the allowed many- 
particle configurations. A description of our numerical methods for analyzing this problem 
follows in Section 3. Results of the numerical study cover both two- and three-dimensional 
systems, and are presented respectively in Sections 4 and 5. Among other results, we provide 
specific examples of interaction potentials whose classical ground state are configurationally 
degenerate and disordered. Our conclusions and discussion of some remaining issues appear 
in Section 6. 

2. GENERAL RELATIONS 

Suppose that the N point particles reside in a one-, two-, or three-dimensional container 
that is an interval, an x Ly rectangle, or an x Ly x rectangular solid. Further- 
more, suppose that periodic boundary conditions apply. The applicable wave vectors have 
components: 
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(2.1) 
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where j — x, y, z as needed. It is easy to see that the C{k) must obey the following 
properties: 

C(0) = ^NiN - 1) , 
C{k) = C{-k) . 

Furthermore, these collective variables are necessarily confined to the range: 

-lN<C{k)<^N{N-l) (fc^O), (2.3) 

and as they vary over this range they measure the magnitude of density inhomogeneity at 
wave vector k in the A^-particle system. 

Although the number of collective variables is infinite, the A'"-particle system possesses 
only dN configurational degrees of freedom, where d is the Euclidean space dimension {d 
— 1,2,3). Consequently, it is unreasonable to suppose (barring special circumstances) that 
generally all C(A;)'s could be independently controlled. However it is possible, as examples 
in Refs. [5, 6] and in subsequent sections below will illustrate, to specify simultaneously a 
number of the collective variables equal to a significant fraction of dN. In particular, let 
Q be a finite set of the fc's meeting this criterion, and let Co{k) be the target value to 
which C{k) is to be constrained. Of course, each Co(fe) must he in the range specified by 
inequalities (2.3) above. Then consider the following non-negative objective function: 

^r,...rN) = Y,V{k)[C{k)-Co{k)]'' , 

feeQ (2.4) 

V(k) = V(-k) > . 

This continuous and differentiable function of the particle coordinates attains its absolute 
minimum value zero if and only if the C{k) for all fc e Q equal their target values. 

By inserting the definitions (1.3) for the collective variables into Eq. (2.4), one trivially 
has the following 

N N N 

$(ri ...rM) = J2 Yl cos[fc-(r -n)] cos[fe-(r^-rO]-2Co(fe) cos[fe-(r -n]+Co'(fe)} . 

feeQ j<l m<n j<l 

(2.5) 

The right member of this last equation can be resolved into symmetric combinations of four, 
three, and two particle contributions, plus an additive constant: 
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(2.2) 



N N N 

^(n ■ ■ ■ r-jv) = 5Z V4{rj,ri,rm,rn) + ^ V3{rj,ri,rm) + ^V2{rj,ri) +^0 ■ (2.6) 

j<l<m<n j<l<m j<l 

The specific forms of these contributions are as follows: 

VAivj, ri, Vm, Tn) = 2 V{k){cos[k ■ {vj - r/)] cos[A; • (r„ - r„)] 

fcGQ 

+ cos[fe • {rj - Vm)] cos[k ■ {ri - + cos[fc • {vj - r„)] cos[A; • - r^)]} , 

V3(rj, ri, r^) = 2 ^ y(fc){cos[fc • (vj - ri)] cos[fe • (rj - r^)] 
feeQ 

+ cos[fc • {ri - Vj)] cos[fe • {ri - r„)] + cos[fc • (r^ - rj)] cos[fe • (r^ - r^)]} , 

V2{r„ n) = J2 V{k){cos^[k ■ {vj - n)] - 2Co cos[fe ■ (r, - n)]} , 
feeQ 

feeQ 

(2.7) 

If $ is interpreted as a potential energy of interaction for the N point particles, then 
Eqs. (2.6) and (2.7) specify four, three, and two particle interaction potentials operating 
in the system. If classical ground-state configurations for the N particles subject to that 
potential exist for which $ = 0, then those configurations necessarily attain the desired 
target values of the collective variables. 

The remainder of our analysis will be restricted to cases for which the wave vector set Q 
consists of all those k lying within a given distance from the origin: 

Q:0<\k\<K . (2.8) 
In addition, we shall also confine attention to the following specific family of forms for Co{k): 

Co(fe) = -A^/2 + D|fc|" , (2.9) 

where the multiplicative constant D and the exponent a are both non-negative. This choice 
focuses on the behavior of the particle system in the hyperuniform regime. 

Previous work [5, 6] only considered the situation where D — 0, i.e., each of the con- 
strained collective variables was required to be at its minimum value —N/2. In that event. 
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it is possible to utilize a simpler potential energy (objective) function $ that reduces to a 
sum of just two-particle interactions: 

$(ri . . . riv) = Q-^ Vik)C{k) , (2.10) 

fcGQ 

In this expression Q is the system length (L^,), area {L^Ly), or volume {L^LyL^), and 

V{k) = V{-k) > . (2.11) 
It is easy to show [5, 6] that Eq. (2.10) is equivalent to: 

N 

l>(ri . . .rjv) = Yv2{rj,ri) , 

(2.12) 

^2(r'j, ri) = ^ V{k) exp[ik ■ {rj - r/)] . 

fcGQ 

On account of the positivity condition (2.11) on V{k), the absolute minimum of $ will have 
the value: 

min (l>) = -N/2 V V(k) (2.13) 

ri...riv — ^ 

fceQ 

if and only if there exist particle configurations satisfying all of the collective variable con- 
straints. Of course the more elaborate (and more general) formulation defined by Eqs. (2.4)- 
(2.7) is also apphcable to this D — case. 

Suppose that system configurations have been found which succeed in producing the 
desired target values for collective variables over the wave vector set Q. Then these config- 
urations obviously satisfy the same target values over the smaller wave vector set Q' C Q. 
But in view of the fact that Q' entails fewer configurational constraints, one can expect that 
a more inclusive set of A^-particle configurations satisfies those constraints. Indeed that is 
exactly what has been found in our previous one- and two-dimensional studies [5, 6], and 
further two- and three-dimensional examples will be reported in this paper. With respect to 
the objective functions $ and $ whose minimization indicates a solution to the constraint 
problem, the corresponding potential energy interpretations demonstrate the presence of 
configurationally degenerate classical ground states for the N particles, including disordered 
or highly irregular ground-state structures. This ground-state degeneracy phenomenon has 
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been the subject of a recent theoretical study [7] that proceeds in a rather different direction 
than from the collective coordinate perspective presented here. The key idea used by Siito 
to prove a theorem about ground states is that the Fourier transform of the pair potential 
V(k) be nonnegative with compact support, which was first employed in Ref. [5]; see also 
Ref. [6]. 

3. COMPUTATIONAL METHODS 

In all calculations reported below, we assume that the system region Q is constrained 
to a unit square in two dimensions or a unit cube in three dimensions, to which periodic 
boundary conditions arc applied. Special attention is devoted to the choice of the system 
population N for both two and three dimensions. In two dimensions, N has been chosen 
such that all particles could be arranged in the square simulation box in an aligned and 
nearly undeformed version of the triangular lattice. For this purpose, the system occupancy 
N is required to be the product of integers 2pq where the rational number p/q is a close 
approximation to the irrational number 3^/^. The corresponding configuration consists of 2q 
rows of p particles. The same approach was used in selecting the system size in our earlier 
two-dimensional study [6]. In three dimensions, N is selected such that all particles can be 
assembled as an aligned face-centered cubic lattice. Thus, the system occupancy is = 
particles where s is a non-zero integer. 

For any given choice of the independent parameters N, K, D, and a, the majority of our 
numerical studies utilized a random number generator to create an initial configuration of 
the particles inside the simulation space Q. As one would expect, this starting point typically 
produces a large positive value of the objective function of interest. We consider the specific 
case of Eq. (2.4) where V{k) is unity for the set Q under consideration. Introduction of the 
original configuration to the numerical optimization tool MINOP [24, 25] results in a search 
for a particle pattern at the absolute minimum of the objective function. Our earlier two- 
dimensional study [6] involved use of the conjugate gradient method [27] as our numerical 
tool of choice. The greater utility of the MINOP optimization technique for the present 
investigation has been an important advantage that merits brief discussion. 
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FIG. 1: (Color online) Tracking of the potential energy for conjugate gradient and MINOP algo- 
rithms. Top panel: Early steps. Middle panel: Intermediate steps. Bottom panel: Late Steps. It 
should be noted that both algorithms were applied to the same d = 1 minimization problem and 
initial particle configuration. The relevant parameters for the minimization problem are = 200, 
K = Wtt, a = 6, and DK^ = 0.01. 



The three plots in Fig. 1 display the tracked potential energy $ during the course of a 
minimization for both conjugate gradient (CG) and MINOP optimizations. Here, we have 
applied the algorithms to $ minimization for a one-dimensional system with the parameters 
N = 200, K = IOtt, a = 6, and DK^ = 0.01. The final values of the objective function $ are 
typical for several cases examined and clearly indicate that the MINOP strategy is better 
suited to finding a numerically precise solution to the problem than is the conjugate gradient 
approach. This significant disparity can be attributed to details of the multidimensional $ 
"landscape" and to the innate differences between the two algorithms. In CG optimization, 
minimization proceeds in a direction conjugate to the old gradient, i.e., that the change in 
the function gradient be perpendicular to the most recent previous direction of minimization. 
In contrast, the MINOP [24] ("dogleg") strategy is as follows: 

Step Let Aj be a pre-set step bound, Pj be the Cauchy step, and rij be the Newton step, 
i.e., 

Pi = (Ig/ZgfGigJg^ = -Uigi 

where gj, Gj, and Hj are the gradient, the Hessian approximation, and its inverse at 
the ith iteration respectively. 

Step 1 Compute rij and |nj|. If |nj| < Aj, we take Axj = rij. If |nj| > Aj, we compute |tnj|, 
where t is defined in Ref. [24]. If |tnj| < Aj, we take Xj = (Aj/|nj|)nj. 

Step 2 If |tnj| > Aj, we then compute the Cauchy step; and, if |pj| > Aj, we take Ax, = 
-(A./|g,|)g,. 

Step 3 If |Pj| < Aj, we take Axj = (1 — ^j)pj -|-^jinj, where 9i is chosen such that < ^j < 1 
and |Axj| = Aj. 

The MINOP algorithm minimizes a real- valued function of any number of variables based 

on user-provided first derivative and function information. In general, it applies a dogleg 
strategy which uses a gradient direction when one is far, a quasi-Newton direction when one 
is close, and a linear combination of the two when at intermediate distances from a solution. 
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4. RESULTS - TWO DIMENSIONS 



We have performed numerical simulations on a variety of system sizes, random initial con- 
figurations, numbers of constrained vectors {i.e., K cutoff), and choices for the independent 
parameters D and a. Calculations have proceeded to attain high precision for the absolute 
minimum of the objective function of interest, $ or $. Note that some initial configurations 
have not yielded the global minimum of the objective function hypersurface. In such cases, 
it can be inferred that there exist some relative minima along the objective function's hy- 
persurface landscape. However, all cases that are reported below in our analysis involve the 
absolute minimum of the objective function. 

In order to simplify the presentation of our results, we introduce the following parameter 



where M{K) is defined as the number of independently constrained collective coordinates 
and d indicates the system dimension. The parameter x is the ratio of the number of 
constrained degrees of freedom to the total number of degrees of freedom in the investigated 
system and has proved to be a fundamental descriptor in the prior one- and two-dimensional 
studies [5, 6]. 

In the following two subsections, we present the results produced by application of the 
MINOP algorithm to minimization of the objective function $ for two-dimensional particle 
systems. In the first subsection, we demonstrate the abihty to tailor the small-] fc| portion of 

the structure factor associated with point particle systems. In the latter subsection, we study 
the effect of manipulating the structure factor S{k) within the distinctive "wavy crystalline" 
regime reported previously in Ref. [6]. 

4.1. Tailoring the Structure Factor - Two Dimensions 

This section describes our successful attempts to manipulate the structure factor S{k) 
of two-dimensional point particle systems. We subjected the objective function $ with a 
randomly generated initial configuration to the MINOP algorithm in order to evolve the 



M{K) 
dN 
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C(fe)'s toward their target values. Our numerical simulations have involved the use of 168- 
and 418-particle systems. Note that these system sizes are the same as those studied in 
Ref. [6]. In addition, we vary the fc-space range parameter K (27r < K < 407r) and the 
multiplicative parameter D (0.01 < DK"" < 150). We have varied a, a key parameter 
determining the nature of the tailored structure factor. In particular, we have generated 
solutions for the Harrison-Zeldovich [17] spectrum (a = 1) as well as for the a — 4, 6, 
8, and 10 minimization problems. We devoted much of our study to the a — 1 and 6 
minimization problems with less detailed attention to the a = 4, 8, and 10 minimization 
cases. As discussed below, our results demonstrate wide latitude in the capacity to control 
the structure factor for two-dimensional point particle systems. 
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FIG. 2: (Color online) Averaged structure factor plots for the two-dimensional minimization prob- 
lem. The relevant parameters are N = 168, K = 207r, DA'" = 75, x = 0.470238. Top panel: 
Harrison-Zeldovich linear spectrum for small k. Bottom panel: |fc|^ spectrum. Each structure 
factor is averaged over 50 realizations. 



Plots of the system structure factor that have been tailored to fit the linear Harrison- 
Zeldovich and spectra are shown in Fig. 2. Here, the structure factor is derived from 
final configurations via Eq. 1.2 and averaged over 50 independent realizations. For the 
purposes of graphical representation, the structure factor is binned over the reported range 
of fc-space. Note that by construction, the contributions to S{k) below the cutoff K for a 
specific set of parameters (fixed A^, K, D, and a) are identical for each realization. Outside 
the cutoff, the S{k) contributions from each realization deviate irregularly from one another, 
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but their average for |fe| > K shows a weak maximum followed by quick decay to unity as |fc| 
increases. Prom these structure factor plots, it is evidently possible to tailor the structure 
factor to either spectrum. In particular, the respective linear and sextic nature of the plots 
at low k are visually clear. 



FIG. 3: Typical two-dimensional final configurations for the investigated spectra. The relevant 
parameters are N = 168, K = 207r, DK"' = 75, x = 0.470238. Top panel: Harrison- Zeldovich 
linear spectrum. Bottom panel: spectrum. 

Figure 3 displays sample two-dimensional final configurations generated by the $ mini- 
mization problem for the linear Harrison-Zeldovich and spectra. The lack of any visible 
long-range regularity is apparent in the pictured configurations for both cases. However, 
some random point clustering appears to be present in the accompanying configuration for 
the Harrison-Zeldovich spectrum. This differs from the spectrum configuration in which 
an effective repelling particle core appears to be present. Finally, we mention that corre- 
sponding analysis of the and spectra minimization problems yield results 
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consistent with the above observations, providing systematic pattern sequences in both r 
and k spaces as illustrated in Figs. 2 and 3. 

4.2. Effect of |fc|^ Spectrum Imposition on the Wavy Crystalline Regime 

One of the qualitatively distinct regimes that was isolated in our earlier two-dimensional 
collective coordinates study [6] was the "wavy crystaUine" regime. This regime is distin- 
guished by patterns consisting of particle columns that display a meandering displacement 
away from hnearity. For N — 418, we have found the patterns to occur when 0.58 < x < 
0.78 [6]. An example of such a pattern is displayed in the top panel of Fig. 4. The result 
of applying the spectrum minimization to this x interval is our focus in this section. 
As above, we subjected randomly generated particle systems to the MINOP [24] algorithm 
to find the absolute minimum of the objective function Our fe-space parameter K is 
determined by the x range over which the wavy crystalline regime prevails. While focusing 
on the |fe|^ spectrum minimization problem, we have varied the independent coefficient D. 
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FIG. 4: Real space particle patterns for a system of 418 point particles. The C(fc) quantities for 
the wave vectors consistent with x = 0.581339 have been constrained. Top panel: Particle pattern 
in the wavy crystalline regime {D = 0). Middle panel: Particle pattern on which the |fc|^ spectrum 
has been imposed with DK^ = 0.01. Bottom panel: Particle pattern on which the spectrum 
has been imposed with DK^ = 10. 
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The middle and bottom panels of Fig. 4 display 418-particle configurations that result 
from imposing the spectrum on particle systems that at D = lie in the wavy crystalline 
regime. The multiplicative parameter D used in generating the configurations in the middle 
and bottom panels is low {DK^ = 0.01) and high {DK^ = 10), respectively, and the cutoff 
K is near the lower limit of the wavy regime of interest. It is important to note that 
all three particle patterns in Fig. 4 were formed from a common initial configuration. A 
comparison of the three figures reveals distortion and disruption of the meandering nature of 
the reference {D — 0) configuration that intensifies as D increases. Evidently the previously 
documented [6] tendency at D = for increasing x to herd particles towards a crystalline 
arrangement is sabotaged by allowing D to increase. 

5. RESULTS - THREE DIMENSIONS 



In this section, we discuss results for three-dimensional configurations that have been 
subjected to the computational algorithm, MINOP [24] requiring absolute minimization of 
$ and First, we extend our two-dimensional analysis to the tailoring of the structure 
factor of three-dimensional particle systems. Second, wc return to an examination of the 
configurational patterns associated with constrained collective density variables in three 
dimensions. This latter aspect extends our earlier work in one [5] and two [6] dimensions. 
Note that all calculations presented in this section have been carried out to the same high 
precision as for the two-dimensional samples discussed in the preceding Sec. 4. 

5.1. Tailoring the Structure Factor S{k) - Three Dimensions 



Our numerical calculations minimizing $ have involved 256- and 500-particle systems. 
Note that these system sizes conform to the guideline discussed in Sec. 3. We vary the 
fe-space range parameter K (27r < K < 407r), the multiplicative parameter D (0.01 < 
DK°' < 180), and investigate both the a — 1 Harrison-Zeldovich spectrum and the a — 6 
minimization problem. 
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FIG. 5: (Color online) Averaged structure factor plots for the three-dimensional minimization 
problem. The relevant parameters are N = 256, K = lOvr, DK°' = 20, x = 0.334635. Top panel: 
Harrison-Zeldovich spectrum. Bottom panel: | A; |^ spectrum. Each structure factor is averaged over 
6 realizations. 
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FIG. 6: (Color online) Snapshots of three-dimensional final configurations for investigated spectra. 
The relevant parameters are N = 256, K = lOvr, DK°^ = 20, x = 0.334635. Top panel: Harrison- 
Zeldovich spectrum. Bottom panel: jfej^ spectrum. 

As before, structure factors are derived from final configurations averaged over several 
realizations. The linear Harrison- Zeldovich and cases are displayed in Fig. 5. Similarly 
to the two-dimensional study, the tailored structure factor deviates irregularly for \k\ > 
K prior to averaging over a set of realizations. Once again the averaged S{k) exhibits a 
peak just beyond |fc| = now an even stronger feature than in the two-dimensional cases 
shown earlier in Fig. 2. An examination of the plotted structure factors clearly reveals the 
linear and sextic nature of the two curves near the origin. Figure 6 displays representative 
three-dimensional final configurations for the $ minimization problem for both the linear 
Harrison-Zeldovich and spectra. The high degree of disorder is evident in both sample 
configurations. 

A system-size scaling study was performed as part of our investigation. This is relevant 
to the approach to the thermodynamic limit. Specifically, we compared the simulation times 
for two different system sizes of fixed x (appropriately scaled independent parameters A^, 
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and D) for the Harrison-Zeldovich spectrum. We find that doubhng the system size 
increases the computation time by approximately a factor of 10. 

5.2. Minimizing Collective Density Variables C{k) - Three Dimensions 

Our numerical studies for D = concentrated on two system sizes (A^ = 108 and 500), a 
wide range of fe-space constraints i.e. low through high and a variety of initial configura- 
tions (both random and lattice-based). For the wide range of fe-space traversed, the odds of 
hitting relative minima of the total energy function $ was increased for large K. However, 
trajectories converged to the absolute $ minimum for all cases that were included in our 
final analysis. Our findings remain substantially consistent between the two investigated 
system sizes. 



TABLE I: Classification of investigated cases associated with each of the regimes. Note that the 
multiplicative parameter D = Q for this objective function $ minimization problem. 



N 


Disordered Regime 


Crystalline Regime 


108 


X < 0.469136 


X > 0.524691 


500 


X < 0.500666 


X > 0.516666 



In the course of our investigation, we observed two distinct regimes of the final configura- 
tions as X (see Eq. 4.1) was varied: disordered and crystalhne. Table I presents the relevant 
range in x for the two investigated systems and distinguishable regimes. The table indicates 
disorder for low values of x and crystallinity for high values of x- This is expected by analogy 
with the reported results of the one- [5] and two-dimensional [6] articles on collective den- 
sity variables. However, our three-dimensional analysis indicates an abrupt transition from 
disordered to crystalline regimes revealing the lack of an intermediate "wavy crystalline" 
regime as observed in two dimensions. For the 500-particle system and x > 0.516666, we 
have verified that the crystal structure is a face-centered cubic lattice, which for the density 
of the system is consistent with the predictions of Siito [7] . 
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FIG. 7: (Color online) Real space particle patterns for the two distinct regimes for systems of 
500 point particles. Top panel: Particle pattern in the disordered regime. The C(fc) quantities 
consistent with parameter x = 0.171333 have been constrained to their minimum values —N/2. 
Bottom panel: Particle pattern in the crystalline regime. The C{k) quantities consistent with 
parameter x = 0.702666 have been constrained to their minimum values —N/2. 
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FIG. 8: Radial distribution functions for systems of 108 point particles. Top panel: The C{k) 
quantities consistent with parameter x = 0.123457 have been constrained to their minimum values 
—N/2. Bottom panel: The C{k) quantities consistent with parameter x = 0.262346 have been 
constrained to their minimum values —N/2. Each radial distribution function is averaged over 10 
realizations. 

Particle patterns respectively within the disordered (top panel) and crystalline (bottom 
panel) regimes for N — 500 are displayed in Fig. 7. The examples shown are typical for 
the two regimes and are vividly distinct. Further insight into the generated point patterns 
follows from examination of the associated pair correlation functions [18] (see Fig. 8). The 
emergence of an effective repulsive core for increasing values of x is apparent, reminiscent 
of a similar effect observed in Ref. [6] in the two-dimensional case. 
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6. CONCLUSIONS AND DISCUSSION 



The former studies of collective coordinate properties presented in Refs. [5] and [6] were 
restricted to one- and two-dimensional point patterns, and documented the effect of forcing 
sets of the collective variables for |fc| < K to their individual absolute minima. The present 
investigation extends that analysis in two distinct directions, by considering point patterns in 
three dimensions, and by examining the effect of constraining the collective variables around 
the fe-space origin to chosen increments above their absolute minima. Specific assignments 
of the increments that have been considered have the form -D|fe|°, where D > 0, and a — 
1, 4, 6, 8, and 10. Point particle configurations satisfying these various collective variable 
constraints have been obtained to high numerical precision starting from both random and 
from distorted-crystal initial conditions, followed by minimization of appropriate objective 
functions and $, Eqs. 2.4 and 2.10]. Although previous attempts have produced configu- 
rations for which the associated structure factor S{k) oc k°' with a < 4 [17, 28], our analysis 
also yields specific configurations that exhibit S{k) oc /c" with a > 6. In addition to the 
specific cases reported in this paper, our studies reveal that a considerably wider range of 
constrained C{k) patterns can be imposed on many-particle systems. For example, forcing 
C{k) to equal instead of its absolute minimum —N/2, reveals a tendency to produce a 
high degree of disorder in the generated particle systems. 

When D — 0, increasing the fraction x of system degrees of freedom subject to collec- 
tive variable constraints was found in one and two dimensions to drive the point particle 
configurations more and more toward crystalline periodicity. The present extension not 
surprisingly shows that the same qualitative trend applies to three dimensions as well. How- 
ever the "wavy crystalline" regime reported in Ref. [6] to separate the two-dimensional low-x 
disordered regime from the high-^ crystalline regime appears to have no three-dimensional 
analog. As x increases for D — cases in three dimensions, the relevant $ landscape on 
which numerical absolute minimization needs to be carried out develops an increasing den- 
sity of relative minima, thereby inhibiting (but not necessarily preventing) the search for 
qualifying particle configurations. For the few cases considered, increasing D above zero at 
constant x has the effect of disrupting the tendency toward crystalline order. 
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A characteristic feature of the various cases examined, provided that the fraction x is 
below its maximum value to yield valid solutions, is degeneracy of the final particle configu- 
rations. This is obvious when several independent initial configurations for a case considered 
are individually subjected to the same minimization operation, and then produce geomet- 
rically distinct final patterns. The minimized objective functions $ and for D > and 
D — ^ cases respectively, indicate that these configuration sets are degenerate classical 
ground states, including disordered ones, for specific potential energy functions. When $ 
is used, that potential function consists of four, three and two body components, while for 
$ only two-body components arise. This offers a specific constructive method to achieve 
degenerate-ground-state potentials, a subject recently discussed by Siito [7]. 

The present collective variable approach may also supply some insight into the existence 
of potentials whose classical ground states are amorphous, or at least highly irregular. Start 
with a small-x constraint case, satisfied with an irregular point-particle configuration. This 
specific configuration will display specific C(fe) values for the unconstrained range |fc| > K. 
Next, formally expand the constraint radius K to the point where x = 1, regarding the 
already-obtained collective variable values as constraint targets. The irregular configuration 
in hand becomes automatically and trivially a proper solution for this extended problem. 
The resulting potential energy specified by $ or $ then possesses that irregular configuration 
as its unique classical ground state (aside from overall translations). 
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